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Abstract. Many barred galaxies harbor small-scale secondary bars in the center. The evo- 
lution of such double-barred galaxies is still not well understood, partly because of a lack 
of realistic N-body models with which to study them. Here we report the generation of 
such systems in the presence of rotating pseudobulges. We demonstrate with high mass and 
force resolution collisionless N-body simulations that long-lived secondary bars can form 
spontaneously without requiring gas, contrary to previous claims. We find that secondary 
bars rotate faster than primary ones. The rotation is not rigid: the secondary bars pulsate, 
with their amplitude and pattern speed oscillating as they rotate through the primary bars. 
This self-consistent study supports previous work based on orbital analysis in the potential 
of two rigidly rotating bars. We also characterize the density and kinematics of the N-body 
simulations of the double-barred galaxies, compare with observations to achieve a better 
understanding of such galaxies. The pulsating nature of secondary bars may have important 
implications for understanding the central region of double-barred galaxies. 



1. Introduction 

Recent imaging surveys have revealed the fre- 
quent existence of nuclear bar s in a large num- 
ber of barred galaxies, e.g., lErwin & Sparkd 
d2002l) found that double-barred galaxies are 
surprisingly common: at least one quarter of 
their sample of 38 early-type optically-barred 
galaxies harbor small-scale secondary bars. 
They found that a typical secondary bar is 
about 12% the size of its primary counter- 
part. The facts that in ner bars are also seen 
in near-infrared (e.g., [Mulchaev et al.1 Il997t 
lLaine et ai1l2002lh and they are often found in 
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gas-poor SOs indicate that most of them are 
stellar structures. Results from these surveys 
also show inner bars are at a random angle rela- 
tive to the primary bars, implying that they are 
probably dynamically independent structures. 
Shlosman et alJ(ll989h invoked multiple nested 
bars to channel gas inflow into galactic cen- 
ters to feed AGN, in a similar fashion as the 
primary bar drives gas inward. However, re- 
cent work suggests that this mechanism may 
not be as efficient as o riginally hoped (e.g., 
iMacieiewski et al.ll2002l) . 

Simulations offer the best way to under- 
stand double barred systems. However, the de- 
coupled nuclear bars that formed in early simu- 
lations did not last long. For example, the most 
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long-l ived nuclear bar in iFriedli & Martinet 
d 19931) lasted for less than two turns of the 
primary bar, corresponding to about 0.4 Gyr, 
which is far too short to explain the observed 
abundance of nested bars. Furthermore, their 
models usually require substantial amounts of 
gas to form and m aintain these nuclear bars. 
iHeller et all d2007allbh reported that nested bars 
form in a quasi-cosmological setting, but the 
amplitudes of the bars also seem to weaken 
rapidly after most of gas has formed stars 
dHeller et alj2007al) . lPetitpas & Wilson! d2004 
found that 4 out of 10 double-barred galaxies 
contain very little molecular gas in the nuclear 
region. These clues suggest that large amounts 
of molecular gas may not be necessary to main- 
tain central nuclear bars. 

On the side o f orb ital studies, 

iMacieiewski & Sparkd (1 19971 |2000) dis- 
covered a family of loop orbits that may form 
building blocks of long-lived nuclear stellar 
bars (see also IMacieiewski & Athanassoulal 
120071) . Their studies are very important for 
understanding double barred galaxies, but 
their models are not fully self-consistent, since 
nested bars in gene ral cannot rotate rigidl y 
through each other dLouis & Gerhard! 1988). 
So fully self-consistent N-body simulations 
are still needed to check if their main results 
still hold when the non-rigid nature of the bars 
is taken into account. 

Here we demonstrate that long-lived sec- 
ondary bars can form in purely collisionless 
N-body simulations, when a rotating pseu- 
dobulge is introdu c ed in t he model (see also 
Debattista & Shed 120071: IShen & Debattista! 
20091) . The nuclear bars in our work are dis- 
tinctly bars, and do not have a spiral shape. 
We show that the behavior of our models are in 
g ood agreement with the loop o rbit predictions 
of Maciei ewski & S parke (2000). We also ana- 
lyze the photometrical and kinematical proper- 
ties of high resolution models. Our theoretical 
results can also be compared to the observed 
2D kinematics of some double-barred galaxies, 
to achieve a better understanding of the dynam- 
ics of the secondary bars. 



2. Model setup 

Our high-resolution simulation consists of a 
live disk and bulge component. We do not in- 
clude a halo component for simplicity, also be- 
cause secondary bars are very small-scale phe- 
nomena in galactic centers where visible mat- 
ter is dominant. The initial disk has the ex- 
ponential surface density profile and Toomre's 
Q ~ 2.0. The bulge was generate d using 
the method o f iPren dergast & Tomer! dl970l) 
as described in iDebattista & Sellwoodl d2000l) . 
where a distribution function is integrated it- 
eratively in the global potential, until conver- 
gence. The bulge has the mass of Mb = 0.2Md. 
The bulge set up this way is un-rotating, we 
then give the rotation of the bulge by simply 
reversing the negative azimuthal velocities of 
all bulge particles into the same positive val- 
ues. We have checked that systems set up this 
way are in very good virial equilibrium. 

3. Results 

Fig. Q] shows the projected double-barred 
model (at t = 405 when the two bars are 
nearly perpendicular) with an ordinary orien- 
tation: the system is inclined at i = 45° with 
the line of nodes (LON) of ifr nuc = 45° relative 
to the secondary bar major axis. The surface 
density image and contours resemble many ob- 
served double-barred systems, such as NGC 
2950, even though we did not deliberately set 
out to match it. 

Fig. [2] shows radial variations of m = 2 
Fourier amplitude and phase for run D at t = 
400. Fig. [3] shows the ellipticity and position 
angle (PA) profiles of ellipses fitted with IRAF 
for the same data as in Fig.|2](we use log scale 
for radius to be consistent with what observers 
usually adopt). There are four popular meth- 
ods for determining the s emi-major axis ob of 
a bar, as s ummar i zed b y lO'Neill & Du binski 
(2003) and Erwin (2005). For convenience, we 
denote the primary bar as Bl and the secondary 
bar as B2: 

(1) the bar end is measured by extrapolating 
half-way down the slope on the m — 2 am- 
plitude plot (Fig Ek). We find «bi ~ 2.3, 
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Fig. 1. The double-barred model at t = 405 pro- 
jected to i = 45° and i/r mc = 45° with all particles 
shown. The model bears a passing resemblance to 
NGC 2950. 



a B2 ~ 0.4, the B2/B1 bar length ratio is 
about -0.17. 

(2) the bar end is measured when m = 2 phase 
deviates from a constant by 10° (Fig |2j)). 
We find fl B i ~ 2.1, a B2 ~ 0.4, the B2/B1 
bar length ratio is about ~ 0.19. 

(3) the bar end is measured at the 
peak of th e fitted ellipticit y pro- 



files (e.g., Marinova & Jogeel 120071: 



(4) 



iMenendez-Delmestre et al.l 120071) . which 
is shown in Fig [3^. We find a B i ~ 1-7, 
a B2 ~ 0.2, the B2/B1 bar length ratio is 
about ~ 0.12. 

the bar end is measured when the PA of fit- 
ted ellipses deviates from a constant by 10° 
(FigEp). We find a Bl ~ 2.3, a B2 ~ 0.4, the 
B2/B1 bar length ratio is about ~ 0.17. 



Method 1, 2 and 4 yield consistent values 
of the bar lengths and length ratios. We found 
that method 3 tends to give a lower value of bar 
le ngths than the other three m ethods, as shown 
in lO'Neill & Dubinskil(l2003l) . Although these 
methods have some uncertainties in measur- 
ing the bar lengths, the length ratio of the two 
bars is in the range of 0.12 to 0.19 (in par- 
ticular method 1, 2, and 4 give a consistent 



Fig. 2. The radial variations of the m =2 Fourier 
amplitude and phase of all particles for Run D at 
t = 400. 



narrow range of 0.17 to 0.19). This result is 
in good agreement with the typical observed 
length ratio of l ocal S2B systems (median ra- 
tio ~ 0.12, see lErwin & Sparke 2002; lErwinl 
120041: iLisker et alj|2006l) . Note that we expect 
that the length of the secondary cannot be too 
large, otherwise the gravitational torque from 
the primary bar will inevitably twist the sec- 
ondary into alignment if they rotate at different 
pattern speeds. 

Fig. [4] shows the behavior of the az- 
imuthally averaged Q, Q + k/2, and the lo- 
cation of the Lindblad resonances of the 
bars at around t - 400. As shown in 
iDebattista & Shed d2007l) . the pattern speeds 
of the bars, especially that of the secondary, 
vary as they rotate through each other: the sec- 
ondary bar rotates slower than average when 
the two bars are perpendicular, and faster 
when the bars are parallel. The patten speed 
bands shown in Fig. [4] reflect such varia- 
tions. Clearly the pattern speed of the sec- 
ondary bar oscillates much more than that 
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Fig. 3. Ellipticity and position angle as a function 
of semi-major axis of IRAF-fitted ellipses for Run 
D at t = 400. 



of the primary. The primary bar extends 
roughly to its CR radius (~ 2.5), consis- 
tent with the general expectation and is there- 
fore considered a fast bar (e.g., Corsini et all 
120031: iDebattista & Willia ms 2004). The sec- 
ondary bar rotates faster than the primary bar. 
However, the secondary bar is much shorter 
than its shortest Rqr- In addition, even if 
the variation of the pattern speed is taken 
into account, the Rcr of the secondary is not 
very close to the Rmg. of the primary, if we 
use the same naive definition of /?ilr as in 
iPfenniger & Normanl fTggoFl This is inconsis- 
tent with the CR-ILR coupling proposed to be 
a requirement for making secondary bars (e.g., 
Pfenniger & Normadfl990l: iFriedli & M artinet 



1993h . 



1 A cautionary note is that the Rius. read naively 
from Fig. |4]serves just as a visual guide, because the 
/?ilr determined this way is reliable only for weak 
bars, and is questionable fo r our strong bars (e.g., 
Ivan Albada & Sandersll 19821) . 




Fig. 4. Frequencies as a function of radius at around 
t = 400 for run D, calculated based on the az- 
imuthally averaged gravitational attraction. The full- 
drawn line shows the curve of the circular angular 
frequency fl and the dashed curves mark CI ± k/2, 
where k is the epicyclic frequency. The two shaded 
bands show the oscillational ranges of the bar pat- 
tern speeds (the upper band is for the secondary bar 
and the lower one is for the primary). 



We also analyzed the line-of-sight ve- 
locity distribution (LOSVD) by measuring 
the mean velocity v and velocity dispersion 
<t. Departures from a Gaussian distribution 
are parametrized by Gauss-Hermite moments 
dGerhardll993l Ivan derMarel & Franxlll993t 
Bend er et al I ll994l) . The second order term in 
such an expansion is related to the dispersion. 

measures deviations that are asymmetric 
about the mean, while /14 measures the low- 
est order symmetric deviations from Gaussian 
(negative for a 'flat-top' distribution, and posi- 
tive for a more peaked one). 

The most striking feature in Fig. [5] is that 
the twist of the kinematic minor axis (i.e., 
vios = 0) in the secondary bar region is weak 
(see the mean velocity maps in Fig. [5^, Eh). 
The kinematic minor axis is almost perpendic- 
ular to the inclination axis, although there is 
a small but noticeable twisted pinch near the 
kinematic minor axis in the nuclear region. The 
weak central twist is mainly due to the rela- 
tively large velocity dispersion, especially in 
the central region (likewise at t — 20 when only 
the small nuclear bar exists, the stellar twist is 
stronger than at t = 405, but still quite small 
compared to the expected twist in gaseous 
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Fig. 5. Photometrical and kinematic maps of our double-barred model. For each row from left to right are 
the projected surface density, mean velocities, velocity dispersion, h 3 , and h 4 maps, (a): Row 1, the double- 
barred model at t = 405 inclined at 45° with the LON of 45° relative to the primary bar major axis; (b): Row 
2, close-up view of (a). One of the dashed lines represents the line of nodes (45°), while the other dashed is 
the anti-LOS (135°). 



kinematics). On the other hand, the twist of 
the kinematic maj or axis is more promi nent in 
the central region. Moise ev" et al.1 (|2004) found 
the stellar kinematic minor axis hardly twists 
from the PA of the disk in their sample with 
the most reliable kinematics, leading them to 
question whether nuclear photometric isopho- 
tal twists represent bona fide dynamically de- 
coupled secondary bars. We demonstrate that 
an authentic decoupled secondary bar may in- 
deed produce a very weak twist of the kine- 
matic minor axis in the stellar velocity field. So 
a central st ellar velocity map w ithout a strong 
twist as in Mois eev et al.1 (120041) does not nec- 
essarily exclude the existence of a decoupled 
nuclear bar. 

4. Conclusions 

We have analyzed the photometrical and kine- 
matical properties of our high resolution mod- 
els, and contrasted them when with or without 
a secondary bar. This study also compared the 
simulated secondary bars with observations. 

In general the shape of secondary bars in 
our models is reasonable compared to observed 
ones. The length ratio of two bars, determined 
by various methods, is in the range of 0.12 to 
0.19, in goo d agreement with lErwin & Sparkd 
(120021 l2003h . The primary extends roughly 



to its corotation radius, and therefore fits 
the definition of a fast bar (see for example 
lAguerri et al.l 120 03). Although the secondary 
bar rotates more rapidly than the primary, its 
semi-major axis is much shorter than its coro- 
tation radius, even if we take the oscillation of 
the bar patterns speeds into account. We did 
not find evidence of CR-ILR coupling (e.g., 
I Pfenniger & Normanlll990l:lFriedli & Martinet! 
1 1993b in our models. 

We find that the central twist of kinematic 
axes is quite weak even if a secondary bar 
is present, due to the relatively large veloc- 
ity dispersion of stars in the central region. 
This is consistent with the 2D s tellar kinemat- 
ics of secondary bars studied in iMoiseev et al.1 
(2004). We do not find a cr (velocity disper- 
sion) drop for our secondary bar model. It is 
more likely that cr-drops are just the signature 
of newly-formed stars, and it is not necessarily 
a unique feature of double-barred systems. 

The general agreement between our sim- 
ulations and observations of double barred 
galaxies gives us confidence that the simula- 
tions are capturing the same dynamics as in 
nature. This is especially remarkable because 
secondary bars are not merely scaled down ver- 
sions of primary bars, but have distinctly dif- 
ferent kinematic properties. In the absence of 
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self-consistent simulations, earlier orbit-based 
models could not directly confront the chal- 
lenge from observations which found such dif- 
ferences. This demonstrates the advantage of 
finally being able to simulate stellar double- 
barred galaxies, which had been puzzling for 
so long. 
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